The extinction law in high redshift galaxies 
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ABSTRACT 

We estimate the dust extinction laws in two intermediate redshift galaxies. 
The dust in the lens galaxy of LBQS 1009-0252, which has an estimated lens 
redshift of zi ~ 0.88, appears to be similar to that of the SMC with no significant 
feature at 2175A. Only if the lens galaxy is at a redshift of zi ~ 0.3, completely 
inconsistent with the galaxy colors, luminosity or location on the fundamental 
plane, can the data be fit with a normal Galactic extinction curve. The dust in 
the zi = 0.68 lens galaxy for B 0218+357, whose reddened image lies behind a 
molecular cloud, requires a very fiat ultraviolet extinction curve with (formally) 
Rv = 12 ± 2. Both lens systems seem to have unusual extinction curves by 
Galactic standards. 

Subject headings: cosmology: observations — gravitational lensing — dust, ex- 
tinction — galaxies: ISM 



^Based on observations made with the NASA/ESA Hubble Space Telescope. The Space Telescope Science 
Institute is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract 
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1. Introduction 

Precise measurements of extinction curves are almost exclusively limited to the Galaxy, 
the LMC and the SMC, because at greater distances it becomes impossible to obtain the 
photometry or spectroscopy of individual stars needed for accurate extinction law measure- 
ments. Galactic extinction curves are well fitted by parameterized models with Ry ~ 3.1 
(^A = RxE{B-V)), e.g. Savage & Mathis 1979, Fitzpatrick & Massa 1988, GardeUi, Glayton 
& Mathis 1989, hereafter GGM), although lines of sight with dense molecular gas can show 
much higher values (e.g. Jenniskens & Greenberg 1993) for a total range of 2.1 < Ry < 5.8 
(see e.g. Draine 2003 and references therein). In the Galaxy and most of the Large Magel- 
lanic Cloud (LMC), the changes in the extinction law are strongly correlated with changes 
in the width and amplitude of the 2175A feature in the extinction curve. The extinction 
curve of the Small Magellanic Cloud (SMC) and some regions of the LMC is very different 
from the typical Galactic extinction law in having a far weaker or nonexistent 2 175 A feature 
(e.g. Misselt et al. 1999, Gordon et al. 2003). Physically, the extinction law depends on 
the mean size and composition of the dust grains along the line of sight (e.g. Clayton et al. 
2003, Draine Sz Malhotra 1993, Rouleau et al. 1997), so it should not be surprising that it 
varies with the environment. 

There are only fragmentary data on the extinction curves in other local galaxies. A range 
of Rv are found in M31, and the variations may be correlated with the local metallicity 
(Hodge & Kennicutt 1982, lye & Richter 1985). There is evidence that the extinction 
curves of early- type galaxies are steeper functions of A~^, but there are no cleanly measured 
extinction curves (Warren-Smith & Berry 1983, Brosch & Loinger 1991, Goudfrooij et al. 
1994). At least in the optical (I-band through B-band), Riess et al. (1996) used Type la 
supernovae to show that the extinction curves of nearby galaxies were consistent with a mean 
optical extinction curve having Ry = 2.6 ± 0.3. In short, aside from one spiral galaxy and 
two irregulars (the Galaxy, the LMC and the SMC) we have few quantitative measurements 
of dust properties. 

The galaxies, that have been studied so far, are not a representative sample of galaxies 
or environments. Moreover, all the physics governing dust properties (metallicity, star forma- 
tion and evolution rates, radiation backgrounds) evolve strongly with redshift, so we would 
expect the properties of the dust to evolve with redshift. Evolution in the mean extinction 
law with redshift would be a crucial systematic uncertainty in studies of Type la supernovae 
to constrain the cosmological model (see, e.g., Perlmutter et al. 1999), since extinction mod- 
ifies the apparent stretch of the light curves (Nugent et al. 2002) . Extinction laws are also 
required for models of galaxy evolution (e.g., in semi- analytic models, Silva et al. 2001, or 
population synthesis models, Gordon et al. 1997), for estimates of star formation rates in 
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individual galaxies (e.g. Pettini et al. 1998; Meurer et al. 1999) or to construct a global 
history of star formation (e.g. Madau et al. 1998; Steidel et al. 1999). Extinction also affects 
the light curves of 7-ray bursts (e.g. Price et al. 2001, Jha et al. 2001), and deriving the 
extinction law from afterglows requires theoretical assumptions about the intrinsic spectrum 
of the burst. With the increasing need for extinction corrections at higher redshifts, it would 
be wise to obtain more quantitative measurements of dust properties at similar redshifts. 
To make any progress, we need a probe of extinction which has the precision of local stellar 
measurements and works at 2; = 1 just as well as at z = 0. 

Gravitational lenses provide a unique tool to study the extinction properties of high red- 
shift galaxies. In most of the ~ 80 known lens galaxies we see 2 or 4 images of a background 
AGN produced by the deflection of light by a foreground lens galaxy. When each image's 
hght traverses the lens galaxy, it is extinguished by the dust at that position. When the 
dust is not uniform, the amount of extinction is different for each image, whose observational 
signature is that the flux ratios of the images depend on wavelength (Nadeau et al. 1991). 
Extinction curves have been estimated for several systems using optical and infrared flux 
ratios (Nadeau et al. 1991, Jaunsen & Hjorth 1997, Motta et al. 2002, Wucknitz et al. 
2003), and Falco et al. (1999) made a general survey of dust properties using the available 
lens photometry. Unfortunately, ground-based observations can study the region around the 
2175A feature only for the highest redshift lenses. The feature is redshifted to wavelengths 
above the atmospheric cutoff (3500 A) for zi > 0.6, and is easily studied only for zi ~ 1. 
Most photometry from Hubble Space Telescope (HST) observations is hmited to the V, I 
and H-bands, since the observations were designed to study the lens and host galaxies rather 
than extinction. In this paper we used near-UV observations with the Hubble Space Tele- 
scope (HST) to study the extinction law of two gravitational lenses near the 2 175 A feature. 
We summarize the observations in §2 and the results in §3. 



2. Observations 

Prom the Palco et al. (1999) lens extinction survey, we selected 3 lenses with signiflcant 
extinction whose redshifted 2175 A feature would lie at longer wavelengths than the Lyman 
limit of the source quasar (i.e. 2175(1 + zi) A> 912(1 + Zg) A). We used a total of 10 orbits 
to observe B 0218+357 (6 orbits), LBQS 1009-0252 (2 orbits) and Q 2337+0305 (2 orbits). 

Table 1 shows a log of our WFPC2 observations. Each image was composed of dithered 
but not CR-split sub-exposures. The original observing request was reduced, and our sub- 
sequent decision to obtain data for all the proposed targets forced us to use only 2 sub- 
exposures for most of the observations. Unfortunately, for Q 2237+0305, we found the short 
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integrations insufficient to produce a useful extinction measurement and we include only 

the photometric results for this system. The sub-exposures were combined using standard 
methods (e.g. as in Lehar et al. 2000) but with additional manual masking to control the 
cosmic rays. Table 2 presents the magnitude measurements for all three systems. 



If mo (A) is the intrinsic QSO spectrum expressed as magnitudes at observed wavelength 
A, then the spectrum of lensed image i, mi(A), is 



where Mj and Ei = E{B — V) are the magnification and extinction of image i, and i?(A/(l + 
zi)) is the extinction curve redshifted to the lens redshift zi. By measuring the magnitude 
differences as a function of wavelength for each image pair (e.g., A and B) 



we can measure the relative magnifications AM — Mb — Ma, extinction differences AE{B — 
V) — Eb{B — V) — Ea{B — V), and the mean extinction curve R(V) without needing to 
know the intrinsic spectrum too(A). We assume that the shape of the source spectrum does 
not vary with time, that there is no wavelength dependence Mi{X) to the magnification due 
to microlensing, and that extinction curve is the same for all images. We will discuss these 
assumptions further below. 

We used a standard statistic to fit the model to the measurements and to determine 
AM, AE, Ry and their uncertainties. We used cither the Cardelli et al. (1989) parameter- 
ized models for the Galactic extinction curve or the Fitzpatrick & Massa (1990) model with 
its parameters set to the values found by Gordon et al. (2003) for the average extinction in 
the SMC. The Galactic models have a strong 2175A feature while the SMC models do not. 
We also attempted to determine the dust redshift za by varying the lens redshift zi in our 
fits (Jean & Surdej 1998, Falco et al. 1999). If the wavelength dependence of the flux ratios 
is due to extinction, then we should flnd that z^ is consistent with zi. Table 3 shows the 
results for these parameters. 

Prom the colors of the lens galaxy in LBQS 1009-0252 or its location on the fundamental 
plane, we estimated that the lens redshift is zi — zpp — 0.881°;°! (Kochanek et al. 2000). 
To date, these estimated redshifts have always been confirmed by subsequent spectroscopic 



3. Analysis and Discussion 
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measurements. An example at similar redshift and with similar photometric data is HE1104- 
1805, where the Lidman ct al. (2000) spectroscopic redshift of zi = 0.73 exactly matched the 
prior prediction of zpp = 0.73±0.04 from the fundamental plane. Fig. 1 shows the magnitude 
differences as a function of the observed wavelength, and there is no sign of a 2175A feature 
at the zpp = 0.88 redshifted wavelength of A""*^ ~ 2.45/im~^. For zi = Zpp the best fit 
with a Galactic extinction curve has = 32.7. In these fits we included the H-band fiux 
ratio from Lehar et al. (2000), although the result changes httle if it is excluded. Galactic 
dust is permitted if the redshift estimate is wrong, as we find a perfect fit (x^ = 0.14) with 
a Galactic extinction law for Zd — 0.31 ± 0.09. With this low redshift, the 2175A feature 
lies outside the range of our data (see Fig. 1). This is very unlikely given the colors and 
structural properties of the lens galaxy because it is almost impossible for a lower redshift 
galaxy to mimic the colors and structure of a higher redshift early-type galaxy (see Lehar 
et al. 2000 and Kochanek et al. 2000). 

The data are, however, very well fitted (x^ = 0.65) by an SMC extinction curve red- 
shifted to the expected zi = zpp = 0.88 (see Fig. 1). The SMC extinction curves lack the 
features needed to estimate a dust redshift, and we can find a good fit for almost any lens 
redshift (0 < 2;^ < 2 at 2cr). Chromatic microlensing can also produce wavelength-dependent 
fiux ratios (see e.g. Yonehara et al. 1999), but the lack of significant changes in the V and 
I flux ratios between 1999.01 and 1999.11 argues against a significant contribution from mi- 
crolensing. In order to obtain such a rapid change in the flux ratio with wavelength using 
microlensing, the source would have to lie in a highly magnifled region of the microlensing 
magniflcation pattern where the time scales would be relatively short. This is easily checked 
by comparing the emission line and continuum flux ratios in a spectrum, since extinction 
has the same effect on both spectral components while microlensing primarily affects the 
continuum ffuxes (e.g. Wucknitz et al. 2003). 

The second lens, B 0218+357, has a known spectroscopic redshift of zi — 0.6847 (Browne 
et al. 1993, Stickel & Kuhr 1993). We combined the 6 WFPC2 ffux ratios with the radio 
ffux ratio of ttlb — tra = 1-40 ± 0.03 (Biggs et al. 1999). The radio ffux ratio provides a 
direct constraint on the true magnification ratio AM because it is unaffected by extinction. 
The magnitude differences (see Fig. 2) show a feature at the redshifted wavelength of the 
2175A feature, but the overall pattern does not correspond to a standard Ry — 3.1 Galactic 
extinction curve. The data are well fitted by the CardeUi et al. (1989) model (x^ = 0.72) for 
Rv — 12 ±2. This is an extrapolation well past the upper values of Ry — 5 actually observed 
in the Galaxy (see e.g. Clayton et al. 2003 and references therein), but it maintains the 
overall structure of the high Ry extinction curves observed in the Galaxy. In this case, the 
dust redshift of Zd = 0.70 ± 0.06 is consistent with the spectroscopic redshift. Moreover, for 
B 0218-1-357 we did not include the Lehar et al. (2000) photometry, yet the model passes 
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almost exactly through the H-band point lying between the radio and optical data. This is 
consistent with local observations that the near-IR extinction curve is universal (e.g. Martin 
& Whittet 1990). 

B 0218+357 is known to lie behind a molecular cloud from the existence of molecular 
absorption lines in the radio continuum (see e.g. Combes & Wiklind 1998). Given Ay ~ 

4 mag, an atomic Hydrogen column density of N{Hj) ~ 10^^ cm~^ (Carilli, Rupcn & Yanny 
1993) and molecular Hydrogen column density of between N{H2) ~ 2 x 10^^ cm^^ ( Wiklind 
& Combes 1995; Gerin et al. 1997) and 5 x 10^^ cm'^ (Wiklind & Combes 1995; Combes 
& Wiklind 1997; see also Mentcn & Rcid 1996), wc find a gas to dust ratio of Nn/Ay ~ 
(1 — 25) X 10^^ cm~^ mag^^. The uncertainties are driven by the large range of the molecular 
column density measurements. Gas-to-dust ratios in the Galaxy have an average value of 
Nh/Av ~ 1.87 X 10^1 cm-2 mag"^ (for Ry = 3.1, Bohhn et al. 1978) with a remarkable 
small scatter ~ 30%. At higher Ry the value of Nn/Ay appears to be greater than the 
average, but the highest value measured, corresponding to p Ophiuchus {Ry = 4.2), is only 
twice the average (see e.g. Kim & Martin 1996). Further study will be needed to understand 
a potential relation between the extreme value of Ry ~ 12 with the huge measured Nn/Ay 
ratio for B 0218+357. 

We explored whether adding dust with a different extinction law in front of the bluer 
image, B 0218+357B, would allow solutions with less extreme extinction curves. These 
models arc under constrained, so we computed models for fixed ratios of the extinction 
between the two images {Eb{B-V) = 0.25Ea{B-V), 0.50Ea{B-V) and 0.75Ea{B -V)). 
The results for a particular test, using a standard Galactic Ry = 3.1 extinction law for the B 
image, are shown in Fig. 3. As we increase the extinction for image B, we find fits with lower 
values of Ry for image A. This behavior is generic over a wide range of choices (2 < Ry < 5) 
for the B image extinction curve. In the final analysis, however, this approach does not seem 
to be a plausible solution because models with enough extinction in the B image to allow 
significantly smaller values of Ry also imply an intrinsic source spectrum with an inverted 
21 75 A feature. This is not physically plausible, so it is unlikely that the B image can suffer 
significant extinction. These conclusions do not change if we assume an SMC-like extinction 
law for image B rather than a Galactic law. 

In summary, both the LBQS 1009-0252 and B 0218+357 lens systems seem to have 
unusual extinction curves. This is not true of all lenses. Falco et al. (1999) found several 
systems with extinction curves reasonably consistent with standard Galactic models, and 
Motta et al. (2002) conclude that the lens SBS 0909+532 is fitted well by a standard Ry = 
2.1 ± 0.9 extinction law. There is one caveat on our results, arising from the fact that we 
ignored the possibility of wavelength dependent magnifications due to microlensing from 
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the stars in the lens galaxy (see e.g. Yonehara et al. 1999). Because the deviations from 
normal extinction curves are large, both systems would have to have large microlensing 
magnifications and correspondingly small source sizes. At least for LBQS 1009-0252 this 
can be checked by studying the flux ratios of the images in the emission lines and the 
continuum of the quasars, where they should be the same if dust is responsible for the 
wavelength dependence of the flux ratios and very different if microlensing is responsible. 
Unfortunately, the source in B 0218+357 is a BL Lac object and lacks strong optical emission 
lines. The flux ratios should also show time variability on time scales of 1-10 years if the 
cause is microlensing, and the lack of significant changes between our present results and 
Lehar et al. (2000) argues against microlensing. 

Our results confirm the unique value of gravitational lenses for studying extinction laws 
at cosmological distances. Since the number of known lenses is increasing rapidly, the only 
barrier to further studies of extinction laws is the need for near-UV HST observations to 
study the region near the 2175A feature. 
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Fig. 1. — The magnitude difference as a function of observed wavelength for LBQS 1009- 
0252. The filled squares are the measurements from these observations and the open squares 
are the earlier measurements from Lehar et al. (2000). The solid hne shows the best fit 
Galactic (CCM) extinction curve, where a redshift of Zd ~ 0.3 is required to avoid predicting 
a 2175A feature in the data (the dotted line shows a CCM extinction curve for a fixed 
Zi = 0.88). The dashed line shows the best fit SMC extinction law at the best estimate for the 
lens redshift {zi = zpp = 0.88) from the color, luminosity and location on the fundamental 
plane of the lens galaxy. The infrared data from Lehar et al. (2000) were included in the 
fits, although the result changes little if we use only the data from the current observations. 
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Fig. 2. — The magnitude difference as a function of the rest frame wavelength for 
B 0218+357. The filled squares are the measurements from these observations, the open 
squares are the earlier measurements from Lehar et al. (2000), and the triangle is the radio 
flux ratio (Biggs et al. 1999). The solid line shows the best fit Galactic (CCM) extinction 
curve. The Lehar et al. (1999) photometry (open squares) was not included in the fit, yet 
the model passes almost exactly through the H-band (Aq ^ = 1.06/xm~^) point lying between 
the radio and optical data. 
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Fig. 3. — Examples of models for B 0218+357 that include extinction of the B image with 
a different {Ry = 3.1) extinction curve than is used for the A image. The contours are 
shown from Xmin + 1 to xLin + ^ assuming that Eb{B — V) was fixed to be 0.75, 0.50 and 
0.25 times Ea{B - V). 



Table 1. Log of WFPC2 Observations 



TARGET DATE-OBS FILTER EXP POS. ANGLE 

(yyyy-mm-dd) (sec) (deg E of N) 



B 0218+357 


2000- 


■03- 


■08 


F300W 


3x700 


100.908 


B 0218+357 


2000- 


■03- 


■08 


F380W 


6x1300 


100.908 


B 0218+357 

J—/ \J ^ W 1 '.J '.J 1 


2000- 


•03- 


-08 


F439W 


3x700 


100 908 


B 0218+357 


2000- 


•03- 


-08 


F555W 


2x2300 


100.908 


B 0218+357 


2000- 


03- 


-08 


F675W 


2x180 


100.908 


B 0218+357 


2000- 


■03- 


-08 


F814W 


2x120 


100.908 


LBQS 1009-0252 


1999- 


•11- 


-05 


F336W 


2x500 


-29.5607 


LBQS 1009-0252 


1999- 


11- 


■05 


F380W 


2x350 


-29.5607 


LBQS 1009-0252 


1999- 


■11- 


■05 


F439W 


2x180 


-29.5607 


LBQS 1009-0252 


1999- 


■11- 


■05 


F555W 


2x60 


-29.5607 


LBQS 1009-0252 


1999- 


•11- 


-05 


F675W 


2x60 


-29.5607 


LBQS 1009-0252 


1999- 


•11- 


-05 


F814W 


2x60 


-29.5607 


Q 2337+0305 


1999- 


•10- 


-20 


F218W 


2x180 


123.703 


Q 2337+0305 


1999- 


•10- 


-20 


F255W 


2x160 


123.703 


Q 2337+0305 


1999- 


•10- 


-20 


F300W 


2x120 


123.703 


Q 2337+0305 


1999- 


■10- 


■20 


F336W 


2x60 


123.703 


Q 2337+0305 


1999- 


■10- 


■20 


F439W 


2x120 


123.703 


Q 2337+0305 


1999- 


■10- 


■20 


F555W 


2x60 


123.703 


Q 2337+0305 


1999- 


■10- 


■20 


F675W 


2x60 


123.703 


Q 2337+0305 


1999- 


■10- 


■20 


F814W 


2x60 


123.703 



Table 2. Photometry 



Lens 



Image 



F218W 



F255W 



F300W 



F336W 



F380W 



F439W 



F555W 



F675W 



F814W 



Q 2237+0305 


A 


17.47±0.13 


16.56±0.08 


16.25±0.19 


15.97±0.23 


17.00±0.10 


16.95±0.15 




B 


18.65±0.33 


17.96±0.22 


17.40±0.22 


17.08±0.45 


17.97±0.44 


17.81±0.34 




C 


18.96±0.32 


18.24±0.27 


17.72±0.29 


17.48±0.20 


18.46±0.32 


18.32±0.42 




D 


19.38±0.14 


18.70±0.21 


18.10±0.07 


17.64±0.27 


18.85±0.08 


18.57±0.14 



LBQS 1009-0252 



18.28±0.05 18.78±0.09 18.77±0.02 18.41±0.07 18.15±0.01 17.93±0.01 

(18.47±0.04) (17.95±0.04) (16.63±C 

22.28±0.17 22.28±0.01 21.94±0.08 21.24±0.09 20.57±0.02 20.11±0.04 

(21.05±0.07) (20.00±0.06) (18.20±C 



B 0218+357 


A 


24.40+0.50 


24.46+0.30 


24.34+0.30 


23.57+0.14 
(23.28+0.20) 


22.78+0.16 


22.01+0.10 
(21.83+0.11) 


(17.52±f 




B 


23.53+0.30 


22.83+0.20 


22.59+0.20 


21.09+0.03 
(21.11+0.07) 


20.15+0.01 


19.42+0.02 
(19.39+0.05) 


(16.94+C 



Note. — 



The magnitudes in parenthesis come from earlier HST observations from Lehar et al. 2000. 
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Table 3. Dust Properties using a CCM Extinction Curve 



Lens 


R(V) 


AE{B - V) 


^dust 




AM 




LBQS1009-0252 


2.5±0.3 


0.41±0.04 


0.31±0.09 


(0.88+°:?^)* 


1.3±0.1 


0.14 


B0218+357 


12±2 


0.30±0.04 


0.70±0.06 


0.6847 


1.4±0.3 


0.72 



Note. — * The spectroscopic redshift for the lens LBQS 1009-0252 is still unknown. 
We use the Fundamental Plane redshift estimation from Kochanek et al. (2000). 



